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Abstract 

Observation of the fine structures (anisotropies, polarization, spectral distortions) 
of the Cosmic Microwave Background (CMB) is hampered by instabilities, 1// noise 
and asymmetries of the radiometers used to carry on the measurements. Addition 
of modulation and synchronous detection allows to increase the overall stability 
and the noise rejection of the radiometers used for CMB studies. In this paper we 
discuss the advantages this technique has when we try to detect CMB polarization. 
The behaviour of a two channel correlation receiver to which phase modulation 
and synchronous detection have been added is examined. Practical formulae for 
evaluating the improvements are presented. 
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1 Introduction 



The fine structures (spatial anisotropies, spectral distortions, residual polar- 
ization) of the Cosmic Microwave Background (CMB), relic of the Big Bang, 
are among the most powerfull tools available for probing the evolution of 
the Universe (for a general discussion see for instance [Partridge 1995] and 
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[Staggs et al. 2001]). Their detection can in fact be used to go back at least 
up to redshisft Z ~ 10 6 — 10 7 , when the Universe was extremely young and the 
matter condensations and objects we observe today not yet formed. So far how- 
ever only spatial anisotropies have been discovered ([Smoot 1992]) and are cur- 
rently studied (e.g. [De Bernardis 2000], [Hanany 2000]). Spectral distortions 
and polarization escaped so far detection and only upper limits to their am- 
plitude have been obtained (see for instance [Sironi 2000], [Sironi et al. 2001], 
[Staggs et al. 1999]). In fact the expected signals are extremely faint when 
compared with spurious effects produced by small instabilities of the receiver, 
1/f noise, pick up of tiny fractions of undesired signals, deviations of the 
system components from their ideal behaviour etc. Therefore many radiome- 
ters succesfully used for standard radioastronomical observations become use- 
less when applied to search for the CMB fine structures: ad hoc systems are 
necessary. For polarization studies correlation receivers are usually preferred 
(e.g. [Sironi et al. 1998], [MAP 2001], [SPOrt 2001], [Torbet et al. 1999], 
[Hedman et al. 2001], [Keating et al. 2001]) because they make possible si- 
multaneous measurements of pair of Stokes parameteres, (U and Q or U and 
V), and in principle allow to detect signals of few fiK. Unfortunately such 
a sensitivity is barely sufficient because the expected amplitude of the CMB 
polarized component is probably below the \iK level. Therefore the receivers 
so far used for studies of the CMB polarization have to be improved. In the 
following we discuss the limits of a standard two channel correlation receiver 
and the improvements in noise rejection and offset cancellation one obtains 
adding phase modulation (at the system front end) and synchronous detection 
(at the back end). 



2 Radiometers and noise 



A radiometer (see figure 1) is a chain of linear devices plus a square law 
detector which amplify and convert the signal s(t), collected by the antenna, 
and the noise n(t), produced by the system components, into DC signals V s (t) 
and V n (t) proportional to the power content of s and n 

V(t) = [\s(t)\ 2 + \n(t)\ 2 ] G = [v s (t) +v n (t)} G = V s (t) + V n (t) (1) 

Here the power gain G includes the detector responsivity, (conversion factor 
between power and output voltage (current)), s and n are electric (magnetic) 
fields with zero mean values, while Vi and Vj are voltages, proportional to the 
power content of the signals, whose mean value is greater than zero: all of 
them fluctuate and behave as noise [Van der Ziel 1954] . 
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Fig. 1. Block diagram of a radiometer - H = radiation collector, HFA = predetection 
amplifier, SWD = Power detector, LFA = postdetection amplifier, r = integration 
time, s = signal, n = noise, s 2 = signal power, n 2 = noise power, G = power gain 



Let's call x(t), E{x(t)} and \i x = E{x(t)} one of these signals, its ex- 
pectation value and its mean value respectively. The signal variance is (e.g. 
[Van der Ziel 1954], [Rohlfs 86]): 

oo 

ol = E{x 2 } - E 2 {x} = J w x (u) dv - (2) 





where w x {y) = \a(u)\ 2 is the signal power spectrum and 



+oo 



a(u) 



x(t)e 



-2-KVt 



dt 



(3) 



the Fourier transform of x(t). In practice we can write 



oo ^max V T 

°i + A = Jw x (v)dv ~ J w x {v)dv= j w x {v)dv (4) 

l/T 



where v min ~ l/T and v max ~ 1/r are the minimum and maximum frequencies 
of the signal fluctuations accepted by the system, r the sample collecting 
time, N the number of samples and T = Nr the total observing time. To 
guarantee that the samples are statistically independent r must be longer 
than the system time constant r' (r > 3 r'). 



Two classes of noise n are considered here: 
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a) white noise (also called random, or gaussian, or steady state). The power 
spectrum is frequency independent: 

w wn (u) dv = A wn dv (5) 



therefore 



2 2 11 Atm 

fion + fJ> x — [A wn (v max — V m in)\ = A wn ( — ) >N^oo (6) 

r 1 r 



No matter which value /i x assumes (in many situations \i x = 0), in case of 
reasonable statistics the noise variance approaches quickly a constant value. 
In this case the rms fluctuations of the mean value decrease as iV and T in- 
crease. Therefore when white noise is dominant one can improve the quality of 
the data collected by a radiometer extending the observing time or increasing 
the number of independent data samples collected. 

b) 1/ f noise. The power spectrum is a power law 

w 1/f (v)dv = ^g- du (7) 



with spectral index a ~ 1. 

-o^i A 1/f ln^ = A 1/f In- = A v/ In TV (8) 



It follows that, when 1// noise is important (a > 1), increasing the observing 
time or the number of samples collected does not help. In fact as T increases 
a growing fraction of low frequency noise is added to the system output whose 
level starts to fluctuates at very low frequencies. This effect cannot be cured 
improving the stability of the system temperature or the performance of the 
power supply. 



3 Application to polarimetry 



3. 1 General layout of a correlation polarimeter 



A common configuration used in radioastronomy for polarimetry is the two 
channel correlation receiver shown in figure 2. Fed by a corrugated horn, an 
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Fig. 2. Two channel Correlation Receiver - H = Horn, IP = Iris Polarizer, OMT 
= Orthomode transducer, LNA = predetection amplifiers, PD = Phase Discrimi- 
nator, DA = Differential amplifiers, r = post detection amplifiers and integrators, 
E, Ei, £2, A, B predetection signals, U\, U2, Us, U4, Si, S2, Oi, O2 postdetection sig- 
nals (see text) 

orthomode transducer (OMT) splits the high frequency signal collected by 
the antenna into orthogonally polarized components with a well defined phase 
difference: = if the components are linearly polarized and = it/ 2 if the 
components are circularly polarized. Inserting or removing an iris polarizer 
between horn and OMT we can set <fi = it/ 2 or <fi = 0. The signals available at 
the OMT outputs are then amplified by separate receivers and finally injected 
into the Phase Discriminator (PD), a network of four Hybrid circuits and 
square law detectors which combines phases and amplitudes of the incoming 
signals ([Sironi et al. 1998], [Peverini et al. 2001]). 

To outline the behaviour of the two channel correlation polarimeter we go to 
the frequency domain. If A(v) and B{y) are the monochromatic signals which 
arrive at the inputs of the Phase Discriminator (PD), and 7 (t) = <p{t) + 
their phase difference, (0 is a constant phase difference which accounts for 
differences between the electrical lengths of the receivers), the PD outputs 
are: 

"E/il \h[A 2 + B 2 + 2AB cos( 7 )] 
U 2 h 2 [A 2 + B 2 - 2AB cos( 7 )] 

U 3 ~ h 3 [A 2 + B 2 + 2AB sin( 7 )] 
U 4 \ [h 4 [A 2 + B 2 - 2AB sin( 7 )] 



(9) 



where hi describes the overall gain of receivers and PD components. Differen- 



tial amplification then gives: 









-u 2 




.s 2 _ 




U 3 


-u 4 _ 





Si[A 2 + B 2 } + ABcos( 7 ) ]' 
5 2 [A 2 + 5 2 ]+^ 2 [ A£?sin( 7 ) ]. 



(10) 



where 

















"2(/i 1 + /i 2 )" 


.*2. 




>3 










.2(/l 3 + /l 4 ). 



(11) 



Finally integration of Si and S 2 over a time length r gives the outputs 0\ and 

o 2 . 



For symmetry reasons an ideal receiver should have hi = h 2 ,hs = conse- 
quently Si should be zero, the constant (7 independent) terms should vanish 
and Si and S 2 should be sinusoidal functions of 7 with zero average value. 
Because it is well known (see for instance [Kraus 1966]) that < A5cos(7) > 
and < AB sin(7) > are linear combinations of the Stokes Parameters, we can 
write: 



or 







-<s x >- 




[o 2 _ 




<s 2 >_ 




or 




" < Si > " 




o 2 _ 




.<s 2 >_ 





K(Qcos(6) -Usm(0)) 
K(Qsm(6) + Ucos(9)) 



K(Qcos(9) -Vsm(9)) 
K(Q sm(9) +Vcos(8)) 



(0 = tt/2) 



(12) 



(0 = 0) 



(13) 



The Milano Polarimeter ([Sironi et al. 1998]) is an example of two channel 
correlation polarimeter similar to the one we described above. Observations 
made with two prototypes (Mk-1 used in 1994 at Baia Terra Nova (Antarc- 
tica) and Mk-2 used in 1998 at Dome C (Antarctica)) showed however that 
in spite of the stability and sensitivity provided by the correlation technique, 
both prototypes suffered gain variations and the system outputs had offsets 
([Sironi et al. 1997], [Sironi et al. 1998], [Zannoni 2000]). In fact small differ- 
ences between the receiver components give hi ~ h 2 , h 3 — h± instead of 
hi = h 2 and h 3 = /i 4 therefore the constant terms do not vanish completely 
and offsets of the system outputs appear. Easily these offsets are large com- 
pared to the amplitude of the sinusoidal terms to be measured. Even worse 
1/ / noise and variations of the offset level caused by gain instabilities, mimic 
signals produced by polarized sources. 

To cure these effects receivers can be enclosed in a (modulation - synchronous 
detection) loop, a technique widely used by radioastronomers (see for instance 
the classical Dicke Receiver ([Kraus 1966])). We can modulate the power signal 
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SQG 




Fig. 3. Modulated Correlation Polarimeter - H = Horn, IP = Iris Polarizer, OMT 
= Orthomode transducer, PhS = Phase shifter, CRX = Correlator, SQG = Square 
Wave Generator, . synchronous detector, r = time integrator 

or the wave signal, in amplitude or, when applicable, in phase. Modulation 
is used to mark the signals we want to detect and produces a shift of the 
average receiver output. The synchronous detector (also called demodulator) 
picks out only the components of the signals and noise marked by modulation 
and exclude all the remaining components, above all the noise components, 
improving the signal to noise ratio. The demodulator type and the modulation 
technique must be matched. 

For polarimetry, where it is essential to preserve phase and amplitude of the 
wave signal, phase modulation of the wave signal E(t) is a natural choice, 
therefore the synchronous detector is a Phase Sensitive Detector (PSD). Fig- 
ure 3 is the block diagram of a Two Channel Correlation Receiver to which 
Phase Modulation and Phase Sensitive Detection have been added. Modulator 
and detector are driven by a periodic signal whose frequency u mod is usually 
between tens and thousands of Hz. 

In the following we analyze the benefits this technique has on the system 
performance. 



3.2 Modulation and Offset Elimination 

Phase modulation of the electric wave E(t) is achieved including into arm 1 of 
the correlation receiver a (0 — 7r), just after the OMT, a phase shifter, driven 
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by a square wave signal of period T mod = \jv mo d (short compared to r), 



r(t) = 



+ 1 nT mod <t <{n + \)T mod 

-1 (n+\)T mod <t< (n + 1) 



(14) 



It multiplies i?i by ±1. An identical phase shifter in channel 2, locked in a 
stable position, equalizes the attenuations in channel 1 and channel 2. 

If E^t) = E ol (t)e lu)t and E 2 (t) = E o2 {t)e i ^ t+<t> \ are the signals available at 
the OMT outputs, the inputs of the Phase Discriminator are 



'A' 








_B_ 




. k 2 E 2 (t) _ 





k x E oX e iu,t r(t) 
k 2 E 02 e^ t+ ^ 



(15) 



where k\ and k 2 account for the system gain between OMT and PD, and 

121 '- r) x k x k 2 \E x E 2 cos{Q + 0)]r(t) 







[5i[|^i| 2 + 


k 2 E 2 \ 


s 2 _ 




M\kxE x \ 2 + 


\k 2 E 2 \ 



(16) 



are the outputs of the differential amplifiers. 

The Phase Sensitive Detector multiplies S\ and S 2 by r'(t), a function similar 
to r(t), and integrates the product for a time r, giving: 



\o 1 - 




-< S 1 x 


o 2 _ 




.< S 2 x 



5i[ Ifci^ip + lfca^al 2 ] <r'(i) > 
<f 2 [ |A; 1 E 1 | 2 + |A; 2 E 2 | 2 ] < r'(t) > 



+ 



+ 



r} x k x k 2 \E x E 2 cos(# + 0)] < r(t) x r'(t) > 
r) 2 k x k 2 \E x E 2 sin(fl + 0)] < r(t) x r'{t) > 



(17) 



Because < r(t) >= and < r'(t) >= now the offset terms vanish even when 
5, ^0(hi^ h l+1 ). 

Three configurations are possible: 

i) system unlocked: r(t) and r'(t) are generated independently therefore 
< r(t) x r'(t) >= 

In this condition (marked by apex ul) 

QUI- 



Qui 



(18) 
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always. 



ii) system locked: PSD and modulator are driven, in phase, by the same func- 
tion r'(t) = r(t) (condition marked by apex /). 



Now < r(t) x r'(t) >=< r 2 (t) >= 1 therefore 



w 2 



r)ikik 2 < E X E 2 cos(# + 0) > 
•q 2 kik 2 < E X E 2 sin(# + 0) > 



(19) 



and one gets the Stokes Parameter (see equations. 12 and 13) 

iii) system locked with a phase difference (time delay At) between the appli- 
cation of r(t) to the modulator and to the PSD (r'(t) = r(t + At), condition 
marked by apex A): 







[o£_ 





ri x k x k 2 < E X E 2 cos(6» + <j>)r(t)r(t + At) > 
r) 2 k x k 2 < E X E 2 sin(6> + 4>)r{t)r{t + At) > 



(20) 



Because T mo d is small compared to r and the time during which the amplitude 
E i and E 02 of E\ and E 2 are expected to vary, we can write 

r]ihk 2 < E X E 2 cos(6» + 0) >< r(t)r(t + At) > " 
r] 2 hk 2 < E X E 2 sin(6» + 0) >< r(t)r(t + At) > . 

Irjik^ < E 1 E 2 cos(9 + 0) >1 , . 

I f] 2 kxk 2 < E X E 2 sin(# + 0) > J 1 j 



where 



1 T 

F = < r(t)r(t + At) > = — J r(t)r(t + At)dt = 



1 - 4 -! — ^ < At < +— 

T v 2 - ~ 2 J 



(22) 



When At 7^ and At ^ ±T/2, |F| < 1 and a reduction of the system 
sensitivity occurs. 



3.3 Synchronous Detection and Noise reduction 



So far we neglected the system noise. Let's now call £j the noise produced 
in channel % and ip r i(f) a random phase. When we add it to the signal Ei 
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equations 15 and 17 become: 



'A' 




_B_ 





fci[^ir(*) +£ 1 e^ 1 W] 
fc 2 [£ 2 e^ +e ) + £ 2 e* 2 ^] 



(23) 



and 



o 2 



-<$i[ I^Exl 2 + |A; 2 £ 2 | 2 + l^l 2 + |A; 2 e 2 | 2 )] < r'(f) > 
.<S 2 [ l^^xl 2 + |A; 2 £ 2 | 2 + l^l 2 + |A; 2 e 2 | 2 )] < r'(t) > 



+ 



+ 



ri x k 1 k 2 \E 1 E 2 cos(9 + <p)} < r(t) x r'(t) > 
r} 2 k x k 2 \E x E 2 sin(# + 4>)} < r(t) x r'(t) > 



+ 



< Nixr 1 + iV* > 

< iV 2 x r' + JV 2 * > 



(24) 



where 



iVi 
N 2 



Si 
5 2 



k\e 2 E 2 + 



V2 



kik 2 Ei(E 2 + e 2 ) 



(25) 



and 



N* 

m 



Si 

s 2 



k\e x E x + 



V2 



k\k 2 e 2 Ei 



(26) 



are noise terms (e\ and e 2 have random phases). As in equation 17 the offset 
terms, (which now include e\ and e|) vanish. 

To appraise the filter action of the synchronous detector we have to evaluate 
the noise variance going back to the time domain, (see equations 2, 3 and 4). 
By integration of equations 24, 25 and 26 over the frequency bandwith of the 
receiver we get the total power measured at the system output: 



'< Vi(*) >" 




<V 2 (t) >_ 





<S 1 (t) + V nA (t)xr(t) + V* 1 > 
< S 2 (t) + V n , 2 (t) x r(t) + Vl 2 > 



(27) 



where Si(t) is the power of the signal, V nyi (t) and V*^(t) are the power of the 
noise associated to JVj and N* respectively. 

We then compute the ratio between the standard deviations of the noise cal- 
culated when the system is locked and unlocked: 



R = 



V nA xr 



+ (Jy* 



+ (Jy* 

n.i 



(28) 
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Assuming the worst conditions (1// noise dominant everywhere) for the power 
spectra of V n ^ and V* i we set wy n i = A^jv and uy* . = A* jv. Moreover 
fci ~ k 2 , £1 ~ e 2 ~ £, E 1 ~ £ 2 ~ £, iVi ~ iV 2 ^"iV, iY7 ~ iV 2 * ~ N*, 
A l ~ A 2 ~ A, A* ~ A* ~ A*, K,i ^ K, 2 ^ K, V;*! ~ V* 2 ~ \/* and 
^ ~ 5 2 ~ (5 << ?? ~ r/i ~ 7] 2 - It follows (N*/N) ~ + e) < 1 and, 

when E is small compared to the system noise, (as usual when we look for the 
fine structures of the CMB), N* « N, A » A* and o v = A In (T/r) >> 
°v* — A* In (T/r). 

To get <7v ni xr first of all we compute the power spectrum wv niXr (i/) of V nji x r. 
Calculations (see Appendix A and [Spiga 2000]) give (here an in the following 
we omit all i indexes) 



4 00 1 



A °° 



7T 2 V (2A; + 1) 2 |^-^| 7r 2 q (2/c + l) 2 \v + z/ fc 



(29) 



Because < x r >= from equation 4 follows: 

1/r 



a V n xr 



J w VnXr (v) dis = J w VnXr (is)dv 



l/T 



1/r 1/r 

4 ^ A r 1 4 ^ A r 1 



dp 



i_ 

7T 2 4r (2fc + i) 



inf^^Um^ fc + 1/ ^ 



(30) 



(l/T < 1/r < z/ fc for each fc). 4 
Finally neglecting k > O 3 terms, we get 



R = — = 



xJ ^o d -l/T \ +l J u mod + l/r 

\Vmod-l/Tj \U mod +l/T 



21nf 



iymax ' ' ^ miri) I ^ mod 
\ In ( Umax I Vmin ) 



(31) 



When u mod - (PSD off) <4 xr -> Aln(T/r) = a 2 ^ because E§° rsiTTF = T 
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Discussion 

Above we got formulae which can be used to evaluate the noise reduction R 
or, equivalently, the improvement of the system sensitivity 1/R one obtains 
adding phase modulation and synchronous detection to a correlation receiver. 
Assuming for instance v m od = 1kHz, r = 5 sec (v max = 0.2Hz) and T = 
60 mm [y m %n = 2.8 10~ 4 Hz), as common in CMB observations, we get R ~ 
5.5 1CT 3 and 1/R ~ 180. In practice the effective reduction can be smaller. In 
fact: 

i) We assumed square wave modulation. It gives maximum efficiency, because 
the modulator reaches almost immediately a well defined status, and keep it 
for almost 50% of the modulation cycle. However the circuits which carries 
out this operation must be carefully studied because spurious signals are eas- 
ily triggered by sharp transitions. For that reason sine wave modulation is 
sometimes preferred. It gives smooth transitions therefore the generation of 
spurious signals can be more easily controlled. However the signal which drives 
the modulator varies sinusoidally therefore the modulator response can vary 
and during an important fraction the modulation cycle be poorly defined. 

ii) No matter which shape is preferred, practical modulating functions, modu- 
lators and detectors are only approximations of the mathematical functions we 
assumed. Deviations of the real components from their model as well as phase 
dependence of the attenuation and impedance of the modulator and detector 
may produce spurious modulations and/or cycle asymmetries. Because when 
present, they give rise to < r >^ and/or < V x r >^ these effects must 
be carefully removed once again through proper design of circuits and devices. 
Residuals which should survive can be cancelled by fine shaping the function 
which drives the synchronous detector (e.g. [Sironi et al. 1990]) 

iii) For technical reasons sometimes the modulator - demodulator loop does not 
include the system front end. For instance in the most recent model (Mk-3) 
of the Milano Polarimeter the front end has been set outside the loop because 
no reliable cryogenic phase shifter was available ([Sironi et al. 2001]). In this 
case to analyze the system we have to split gain and noise in channel 1 in 
two components: koi and e i, gain an d noise of the section which precedes the 
loop, kn and en, gain and noise of the section inside the loop. Now k\ = 
k ik 11 ,e 1 = e 01 + (k 11 /k 1 )e 11 = (a+P)e 1 (a = e 01 /e 1 , (3 = (ku/kjen/ei), 
and equation 23 becomes 

'A' 

_B_ 

By an analysis similar to the one we made before, for the noise terms we get 
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(32) 



now: 



i 

2 . 



N' 



Si 
8 2 . 



[Afcoi^i + k 2 2 e 2 E 2 ] + 









N 2 _ 


< 


N 2 _ 



k n k2£u(E2 + £2) 



(33) 



and 



N'* ± 
N't 



A. 



fclfcll£ll(^l + £ l) + 



+ 



fcifc 2 [Sie 2 + £01 (E 2 + £2) 



N* 



1 + a(l + e/E) > 



No 



(34) 



where, to evaluate the approximated expression, we set S\ ~ e 2 — £ and Ei ~ 
E 2 ~ E. The noise component A^'* unaffected by modulation/demodulation, 
is now larger by a factor ~ 1 + (£oi/ £ i) + { e oi/Ei). To keep it small and 
preserve the efficiency of the modulation/demodulation process, £ i an d k i 
must be small compared to S\ and k\ and, even more important, possibly free 
from 1/f contribution. Viceversa the component of the total noise one can 
control through modulation, N', decreases by (3. 

iv)Last but not least u mod must be far from harmonics of all periodic sig- 
nal signals used into the receiver. Among them the frequency u ac of the AC 
power supply is particularly dangerous: it can in fact induce, through residu- 
als ripples on the DC outputs of the power supply which feed amplifiers and 
active components, modulated signals which are picked up by the synchronous 
detector if v mod is close to harmonics of u ac . 

All the effects we listed above reduce the efficiency of the modulation / syn- 
chronous detection techniques in improving the stability and in rejecting the 
noise of a radiometer. By carefull design of the circuitry which realizes the 
system the degradation of R can however be contained sufficiently to make it 
a second order effect. 

The above results have been obtained assuming a phase modulated correla- 
tion receiver. They can be extended to other radiometer configurations, like 
the classical Dicke Receiver ([Kraus 1966]) or bolometric systems which use 
amplitude modulation. In doing it we must remember that modulation and 
synchronous detection have different effects on the performance of a receiver. 
An important effect of modulation is a shift of the average output of the re- 
ceiver. This shift can be used to bring to zero the average value of the receiver 
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output, making the system insensitive to gain fluctuations and allowing large 
amplifications without saturation. Synchronous detection improve the signal 
to noise ratio, creating a filter which excludes signals not marked by modula- 
tion and cutting the components of the noise at frequency different from and 

below V mod . 

Used in the past for classical radioastronomical observations and in many 
physical experiments (the so called lock in technique) modulation and syn- 
chronous detection are today essential to reach the sensitivities necessary to 
study the fine structures of the CMB. 
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A Appendix A : Power spectrum of the modulated noise (V x r) 



By proper choice of t, the square wave r(t) of period T mod = \jv mod can be 
represented by a Fourier serie containing only cos terms and v mod odd multiples 
v k = (2k + \)v mod : 

a +°° (-l) k 

r(t) = -^A-^ T cos(2vr^) (A.l) 



The noise V can be written as a Fourier integral 

+00 

V(t) = J a(z/)e i(2 ™'+^ (i)) dv (A.2) 



where ip u (t) is a randomly variable phase (here and in the following we will 
omit pedix % which marks the system channel). Therefore 

V(t) x r(t) = / a(u) cos(27rz/t + if> v (t)) dv x - V — — cos(2vrz/ fc t) 

J 7T „ 2k + 1 
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+°° • oo ( i\k 



r 4 00 ( — l) h 

+i / a(v) sin(27ri/t + ip v (t)) dv x - V — ^- cos(2vrz/ fe t) (A. 3) 

J 7T „ 2ft + 1 



Using standard trigonometric formulae we can write 
9 +°° C— \ \ k 

+00 



x { y a(i/) [cos(2vr(z/ + i/ fe )i + + cos(2tt(z/ - v k )t + VUO)] rfz/ + 



+00 



+i J a(u)[sm(2n(u + i/ fc )i + ^„(i)) + sin(27r(z/ - v k )t) + ^„(t)] dz/} (A.4) 



Rearranging terms and omitting for semplicity the time dependence of ip we 
get: 

+00 

J a(z/)[cos(27r(z/ + v k )t + ip v ) + cos(27r(z/ — z/ fc )t + ip v )}dv = 



+00 +00 



y a(]s) cos(2ir(]s + v k )t + tp^dis + J a(u) cos(27r(f — u k )t + tp 1/ )du+ 








+ y a(z/) cos(27r(z/ + z/,t)t + ip^dv + y a(z/) cos(27r(z/ — v k )t + ip u )dv 



+OO +^/c 



y a(z/ - z/ fc ) cos(27rz/t + ip v - Vk )dv + y a(z/ - v k ) cos(27rz/t + i]) v _ Vk )dv+ 



+00 



+ y a(i/ + is k ) cos(27ri4 + ^) V+Vh )dv 



+ y a(z/ + v k ) cos(27rz/t + ip v+Vk )dv (A. 5) 



and 

+00 



y a(i/)[sin(27r(i/ + v k )t + Vv) + sin(27r(z/ - z/ fe )t + ip v )]dv 
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+00 +v k 

j a(u — u k ) sin(27ri4 + ty v - v Adv + J a(v — v k ) sm(2irvt + ip 1/ - 1/k )dv+ 



+v k 
+00 



+ J a(v + v k ) sm(2irist + ip v+Vk )di> 



-Vk 

+ J a(v + v k ) sin(27rz/t + ^ V+Vk )dv 



(A.6) 



Therefore 



V(t) x r(t) 



2 +2? (-l) fc 



+00 

/ a(z/ — u k ) exp[i(27rz/t + ip u _ Uk )]du + 



+ J a(v - v k ) exp[i(2irvt + i> v _ Vk )\dv 



+00 

+ y a(f + i/fc) exp[?(27rz/t + ip u + Uk )}dv+ 



-v k 



+ 



-Vk 

j a(v + v k ) exp[i(2irist + ip v+Vk )]dv 



(A.7) 



Because ip v ±v k is a random function of u, the fc component of the power spec- 
trum of V x r is a sum of power spectra: 



w k (v) 



4 



7T 2 (2fc + l) 2 



\a{—u — v k )\ — 00 < v < —v k 
\a(+u + v k ) | 2 — z/ fe < i/ < +00 



+ 



+ < 



\a(u k — v)\ 2 — 00 < v < +u k 
\a(v — v k )\ 2 + v k < v < +00 



(A.8) 



which, when the noise is 1/ f noise (\a(u)\ 2 = \a(— u)\ 2 = A/\v\) becomes 
4 



w k (u) 



n 2 (2k + iy 



A A 
\v-v k \ \v + v k \ 



(A.9) 
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Because the phases of the Wk terms vary very rapidly in a random way (they 
are calculated at frequencies different for each k) , the total power spectrum is 
obtained adding the incoherent terms 

4°°1 A 4°°1 A 

WV „ (V ) = - - - - _ - + _^ ___ _ (A.10) 
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